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ABSTRACT 


"The  time-differential  perturbed  angular  correlation  (TDPAC)  technique 
has  been  used  to  investigate  the  temperature  dependence  of  the  electric- 
field  gradient  (EFG)  at  Cd  impurities  in  the  semiconducting  compound  In2Se3. 
Between  room  temperature  and  the  melting  point,  four  different  phases  could 
be  observed.  Most  crystal  structures  which  have  been  proposed  for  the  a  and 

8  phases  are  shown  to  be  incorrect.  The  6'  phase  (550-720°C)  has  not  been 
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reported  previously.  The  data  suggests  that  the  3  phase  may  be  a  Peierls 
distortion  of  the  6‘  phase.  The  EFG  in  the  a  phase  is  well  described  by  a 
simple  point  charge  model  suggesting  that  the  electronic  part  of  the  EFG  in 


a  -  In2Se3  is  small . , 

I  .  M  ~  r 

I 

/ 


I.  Introduction 


Like  most  III-IV  compounds  of  the  general  formula  ^2'  *n2^e3  is  a 
semiconductor  with  a  number  of  structural  modifications  and  complicated  low 
temperature  phases.  Although  it  has  been  the  subject  of  several  X-ray  and 
electron  microscopy  studies  (1-7),  the  different  structures  and  their  transi¬ 
tion  temperatures  are  not  all  well  established.  The  time-differential  per¬ 
turbed  angular  correlation  (TDPAC)  technique  is  a  powerful  tool  for  investi¬ 
gating  electric-field  gradients,  magnetic  hyperfine  fields  and  relaxation 
mechanisms  at  the  site  of  radioactive  tracer  nuclei.  Electric-field  gradients 
are  very  sensitive  to  the  local  symmetry  of  the  tracer  nucleus,  and  it  is 
often  possible  to  distinguish  between  structures  which  are  only  slightly 
different  (8).  Our  interest  in  the  compound  In^Se^  initially  arose  from  an 
investigation  of  motional  correlation  times  at  1 1 1  In-1  ^ 1  Cd  probe  nuclei  in 
liquid  Se  using  the  TDPAC  technique.  At  lower  temperature  a  solid  often 
precipitated,  and  it  was  assumed  that  this  solid  was  I n^Seg  (9).  The  necessity 
of  characterizing  the  TDPAC  spectrum  of  this  precipitate  and  the  lack  cf 
definitive  structural  information  about  higher-temperature  phases  of  Ir^Se-j 
prompted  us  to  undertake  this  investigation  of  the  temperature  dependence 
of  the  EFG  at  the  site  of  Cd  impurities  in  ^Se^. 

II.  Experimental  Details,  Data  Analysis  and  Results 

111 

The  radioactive  isotope  In  used  in  these  measurement  has  a  half-life 
of  2.83  days  and  decays  by  electron  capture  to  the  419  keV  state  of  11  "'cd. 

This  excited  state  (half-life  0.12  nsec)  subsequently  decays  via  a  y-y 
cascade  to  the  ground  state  of  ^Cd.  The  relevant  nuclear  parameters  for 
this  system  are  given  in  Table  1. 

The  experimental  setup  for  the  time-differential  perturbed  angular 
correlation  measurement  consists  of  four  Nal(Tl)  scintillation  counters. 
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arranged  in  fixed  positions  symmetrically  in  the  horizontal  plane  around  the 
sample.  Delayed  coincidences  of  the  172-247  keV  cascade  were  detected  for 
interdectector  angles  8  3  90°  and  9  3  180°  by  standard  fast-slow  coincidence 
circuits.  The  resulting  four  spectra  were  stored  in  separate  memory  banks  of 
a  multichannel  analyzer  as  a  function  of  the  time  delay  t  between  the  two 
y-emissions.  After  data  were  collected  for  approximately  24  hours,  the  ratio 

2  fpC7(180°,  t)  •  C4(180°,  tU  ^ 

R(t)  =  3  V-C'3  "(9da,  t)  •  C2(90a,  t)  3  '  1| 


was  computed.  Here  the  c-  (9,  t) ,  i  =  1  to  4,  are  the  time  spectra  collected 
by  the  different  detector  combinations  corrected  for  random  coincidences. 
Representative  R(t)  data  are  shown  in  Fig.  1.  R(t)  is  related  to  the  per¬ 
turbation  function  G22(t)  for  electric  quadrupole  interactions  by  (10) 


R(t)  » 


^22  ^22(T) 

^22  ^22 ( ^ ) 


A22 ( t )  1S  anisotropy  parameter  given  in  Table  1. 

The  ^In  tracer  was  obtained  commercially  as  a  carrier-free  solution 
in  HC1 .  Samples  were  prepared  by  drying  a  small  amount  of  the  radioactive 
solution  on  a  few  pellets  of  Se  and  In  which  were  then  sealed  in  a  2  mm  ID 
quartz  tube  under  vacuum.  The  capsules  were  heated  to  a  temperature  just 
above  the  melting  point  of  Ii^Se^  with  occasional  agitation  for  12  hours. 

Then  they  were  cooled  to  approximately  800°C  and  annealed  for  three  days 
before  data  were  collected.  The  samples  were  deliberately  prepared  with  an 
excess  of  Se  above  the  InjSe^  composition  to  minimize  the  possibility  of 
producing  any  of  the  more  In-rich  phases.  At  temperatures  below  750°C, 
the  Se  solubility  in  Ir^Se-j  and  the  In  solubility  in  Se  is  sufficiently  small 
that  the  excess  Se  should  have  no  effect  on  the  TOPAC  spectra. 


The  observed  R(t)  spectra  were  analyzed  by  performing  least-square 
computer  fits  to  the  following  two-component  perturbation  function  for  an 
electric  quadrupole  interaction  in  polycrystall ine  samples  (11): 

3  2 

G22(t)  =  f  l  S2n  cos  (u>ln  t)  exp(~i  uln  5^  t2) 
n=0 

3  2  ?  ? 

+  (1-f)  l  S2n  cos  U2n  t)  exp(-i  w2n  $2C  t4) 

Here,  f  and  1-f  are  the  fractions  of  the  "'"'"'cd  nuclei  which  are  subject 
to  an  electric  quadrupole  interaction  characterized  by  the  frequencies  'j.  , 
and  w2n  respectively.  The  are  functions  of  the  quadrupole  frequency 
vn  =  eQV  / h  and  the  asymmetry  parameter  n  =  (V  -V  )/V  .  The  V..  are  the 

components  of  the  diagonalized  EFG  tensor  arranged  so  V  <  V  <  V  ,  while 

aa  "  jy  * 

the  5.  describe  the  relative  width  of  the  Gaussian  frequency  distribution 
assumed  for  both  lattice  sites.  This  function  provided  a  good  fit  to  all 
data  below  720°C,  and  no  indication  was  found  that  more  than  two  different 
symmetry  sites  for  In  atoms  exist  in  any  of  the  samples.  The  parameters 
obtained  from  these  fits  are  given  in  Table  2,  and  the  derived  EFG  and  n  values 
are  shown  in  Fig.  2.  Between  room  temperature  and  the  melting  point,  four 
different  phases,  denoted  by  a>  8,  3‘,  and  y  can  be  identified.  Above  840°C, 

R(t)  could  be  fit  by 

G22(t)  =  exp(-U) 

indicating  either  a  time-dependent  interaction  (12)  or  a  very  broad  distribution 
of  small  electric-field  gradients  (11).  The  fit  parameters  X  are  given  in  Table  2. 
In  this  temperature  range  however,  the  composition  of  our  samples  probably 
departs  significantly  from  In2Se2»  so  X  probably  does  not  reflect  the  TDPAC 


spectrum  of  pure  In2Se2- 
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Below  550°C  (a  and  3  phases),  the  data  are  well  fit  with  »  *2  ~ 
implying  two  well-defined  lattice  sites.  The  two  sites  are  equally  populated 
(f=0.5).  The  geometric  factors  S,,n  differ  slightly  from  the  theoretical  values 
for  randomly  oriented  polycrystall ine  samples  (13)  and  were  treated  as  free  fit 
parameters.  Between  550  and  720°C  (3‘  phase)  the  data  were  best  described  by 
a  single  lattice  site  and  a  frequency  distribution  of  approximately  -1%. 

In  the  y  phase  (720-810°C),  ^22 ( ^ )  was  constant,  implying  a  cubic  In  environ¬ 
ment.  The  3  phase  shows  a  strong  supercooling  tendency,  and  the  lifetime  of 
the  metastable  3  phase  above  100°C  was  much  longer  than  the  half-life  of  the 
^In  tracer.  The  ct  phase  samples  were  made  by  quenching  from  800°C  to  room 
temperature  and  annealing  at  180°  to  relieve  strains  as  much  as  possible. 

III.  Discussion 

A.  The  Electric  Fieldgradient 

The  observation  of  well-defined  angular  correlation  curves  for  static 
quadrupole  interactions  indicates  that  the  ^Cd  probe  nuclei  are  in 
equilibrium  with  their  electronic  environment  when  the  y-y  cascade  occurs. 

This  is  expected  to  be  the  case,  since  the  0.12  nsec  half-life  of  the  initial 
state  of  the  cascade  is  long  enough  and  the  electrical  conductivity  of  In2Se3(l4) 
is  large  enough  for  the  electron  cores  of  the  Cd  probe  to  relax  after  the 
electron  capture  decay. 

Electric-field  gradients  at  probe  nuclei  in  solids  are  usually  taken  to 
be  the  sum  of  a  lattice  part  V^,  caused  by  the  ion  cores  of  nearby  atoms, 
and  an  electronic  part  which  arises  from  the  influence  of  covalent  bonds, 
conduction  electrons,  and  unfilled  electronic  shells  of  the  probe  atom  (15). 

Since  In2Se3  is  a  semiconductor  and  Cd  has  no  unfilled  d  or  f  shells,  the 
latter  two  contributions  should  be  negligible.  We  have  made  a  lattice 
sum  calculation  of  for  the  a  phase  using  a  point  charge  model  for  the 


nearest  neighbor  atoms.  The  influence  of  second  and  higher  neighbor  shells 
as  well  as  lattice  distortions  due  to  the  substitution  of  Cd  at  the  In  sites 
has  been  neglected.  We  used  the  structural  data  of  tikforman  et  al  (3)  who 
describe  the  a  phase  as  a  distorted  Wurzite-type  structure  with  two  equally 
populated  In  sites.  No  other  of  the  numerous  structures  which  have  been 
proposed  (5,  7)  has  2  equally-populated  In  sites.  The  charge  state  of  the 
Se  near  neighbors  is  estimated  from  the  Se-In  electronegativity  difference 
to  be  -0.364e  (16,  17).  To  correct  approximately  for  the  polarization  of 
the  Cd  electronic  core,  the  field  gradient  due  to  the  point  charges  is 
multiplied  by  an  antishielding  factor  (1-ym)  =  30  appropriate  to  a  Cd2+ 
ion  (18).  The  results  are  shown  in  Table  3.  If  we  identify  site  1  as  the 
tetragonally  coordinated  site  and  site  2  as  the  pentagonally  coordinated 
one,  the  agreement  between  the  calculated  and  the  experimentally  ob¬ 
served  field  gradients  is  surprisingly  good,  so  any  electronic  contributions 
to  V22  are  apparently  small.  The  asymmetry  parameters  are  in  less  good 
agreement,  but  these  are  very  sensitive  to  small  lattice  distortions  around 
the  Cd  impurity  which  were  ignored  in  this  calculation. 

Similar  calculations  cannot  be  done  at  present  for  the  higher  temperature 
phases,  because  detailed  structural  data  are  not  available.  It  is  known 
however  that  IngSe^  contracts  significantly  in  the  a  ■*  3  transition  (19). 
Qualitatively,  one  expects  a  larger  when  the  atoms  are  closer  together, 
and  this  is  consistent  with  the  larger  V22  found  in  the  3  phase. 

B.  Structure  of  In?Se? 

These  TDPAC  measurements  show,  that  all  ^Cd  tracer  atoms  occupy  one 
of  two  distinguishable  In  sites  in  the  a  and  3  phases  of  InjSe^.  Both 
sites  are  equally  populated  and  apart  from  one  site  in  the  a  phase,  all  have 
less  than  axial  syimietry.  There  have  been  a  number  of  X-ray  and  electron 
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diffraction  investigations  of  In^Se^,  and  several  mutually  inconsistent 
crystal  structures  for  the  a  (3,  5,  7)  and  3  (4,  7)  phases  have  been 
proposed.  For  the  a-phase,  only  the  structure  proposed  by  Likforman  et 
al  (3)  is  even  qualitatively  consistent  with  the  TDPAC  results  as  discussed 
above.  For  the  3-phase,  the  only  structure  which  may  be  qualitatively 
consistent  with  the  TDPAC  results  are  those  proposed  by  Semilitov  (4), 
but  his  structural  investigation  was  too  imprecise  for  any  meaningful 
comparison.  The  3'  phase  has  not  been  previously  reported.  In  order  to 
escape  observation  by  numerous  investigators  the  3  and  3'  phases  probably 
are  structurally  similar,  and  the  6^3'  transition  probably  has  a  relatively 
small  specific  heat  anomaly.  This  observation,  along  with  the  temperature 
dependence  of  the  field  gradients  in  the  3  and  S'  phases  leads  us  to  speculate 
that  3  is  a  Peierls  distortion  of  3'.  8oth,  the  relatively  large  V  distri¬ 
bution  of  the  3'  single  site  and  the  convergence  of  the  two  Vzz  and  r  values 
with  rising  temperature  in  the  3  phase  are  suggestive  of  this  possibility. 

It  should  be  emphasized  that  the  appearance  of  the  3'  phase  is  not  an  arti¬ 
fact  of  the  data  analysis.  The  TDPAC  data  of  this  phase  could  be  fit  as 
well  to  a  two  lattice  site  model  with  only  slightly  different  interaction 
frequencies.  Both  frequencies  fall  within  the  width  of  the  Gaussian  distri¬ 
bution  determined  by  the  one  site  model  and  are  significantly  different  from 
the  values  in  the  3  phase. 

We  find  the  y  phase  to  be  cubic  as  reported  by  Likforman  et  al  (2)  but 
can  obtain  no  further  information  from  TDPAC.  Unfortunately  we  also  can 
deduce  no  structural  properties  from  our  data  above  800°C  because  of  the  Se 
excess  in  our  samples.  In  this  temperature  region  a  monoclinic  5  phase  has 
been  reported  (1).  These  TDPAC  data  are  taken  at  temperature  intervals 
which  are  too  widely  spaced  to  reveal  accurately  the  transition  temperatures. 


y  transition 


The  a  -*■  3  transition  temperature  of  22Q°C  and  the  3'  -*■ 
temperature  of  720°C  reported  by  Likforman  et  al  (2)  are  consistent  with  our 
measurements.  The  3  -  3'  transition  occurs  at  approximately  550°C. 
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Figure  1 


Figure  2 


FIGURE  CAPTIONS 


Representative  R(t)  spectra  at  three  temperatures.  Statistical 
errors  are  shown.  The  solid  curves  are  computer  fits  described 
in  the  text. 

Electric-field  gradients  and  asymmetry  parameters  derived  from 
computer  fits.  The  two  sites  in  the  a  and  3  phase  are  equally 
populated.  The  single  site  in  the  3'  phase  has  a  frequency  dis 
tribution  of  4%.  The  lines  joining  points  have  no  theoretical 
significance . 
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Table  1 :  Summary  of  relevant  nuclear  characteristics  of 
radioactive  probe  nuclide  ^In  -►  ^Cd 


Parent  half  life 

2.83  days 

Energies  of  y-y-cascade 

171-245  keV 

Nuclear  spin  sequence 

l/2+  -  5/2+ 

Angular  correlation  coefficients: 

^22 

-0.180  (2) 

A 

"44 

+0.002  (3) 

Intermediate  nuclear  state: 

half-life 

84  nsec 

spectroscopic  nuclear 

+0.83  (13)  b 

quadrupole  moment 

All  data  taken  from  S.  Raman  and  H.J.  Kim,  Nucl .  Data  B6,  39  (1971)  except  the 
quadrupole  moment  which  is  from  P.  Herzog  et  al ,  Z.  Phys.  A  294,  13  (1980). 
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Table  2:  Quadrupole  frequencies  and  asymmetry  parameters 
derived  from  computer  fits  as  described  in  text. 

Static  quadrupole  interaction: 


T[°C] 

yq[mhz] 

n 

5 

25 

115.8  (4) 

0 

0 

(site  1 ) 

83.5  (3) 

0 

0 

(site  2) 

107 

0.11  (4) 

0 

(site  1) 

0 

0 

(site  2) 

182 

120.4  (3) 

0.11  (4) 

0 

(site  1 ) 

82.3  (3) 

0 

0 

(site  2) 

120* 

143.0  (6) 

0.20  (4) 

0 

(site  1 ) 

98.0  (5) 

0.35  (5) 

0 

(site  2) 

181* 

138.0  (6) 

0.28  (4) 

0 

(site  1 ) 

98.0  (5) 

0.32  (5) 

0 

(site  2) 

274 

C Co  -Ar 

0 

(site  1) 

0 

(site  2) 

341 

0 

(site  1 ) 

0 

(site  2) 

407 

129.3  (3) 

0.26  (3) 

0 

(site  1 ) 

98.7  (3) 

0.30  (3) 

0 

(site  2) 

510 

123.4  (3) 

0.25  (3) 

0 

(site  1 ) 

101.7  (3) 

0.27  (2) 

0 

(site  2) 

584 

BSSif^Si 

0 

3.9  (1)S 

628 

130.3  (4) 

0.08  (4) 

3.4  (1)S 

667 

0 

3.3  (1)S 

716 

1 

0 

3.7  (2)S 

759 

0 

— 

— 

806 
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Table  2  (continued) 


Time  dependent  quadrupole  interaction: 


T[°C1 

xTmhz] 

843 

0.29  (5) 

861 

0.26  (5) 

933 

0.09  (5) 

♦supercooled  sample 


-I  Sa¬ 


lable  3:  Comparison  of  experimental  and 
calculated  EFG  for  the  a  phase 


Site 

V22[10  v/cnr]  | 

i 

1  " 

t 

#1  (experimental ) 

4.17  (65) 

0(T  =  25°C) 

i 

0.1 (T  >  25’C) 

4-fold  coordinated  (calculated) 

4.21 

i 

0.23 

#2  (experimental) 

5.78  (91) 

0 

5-fold  coordinated  (calculated) 

4.67  j 

0.19 
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